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I We previously (G. Coddens, Phys. Rev. 63, 064105 (2001)) derived a theorem about the 

, coherent quasielastic neutron-scattering signal from a d-dimensional lattice of A'' molecules that are 

^ ' undergoing rotational jump diffusion (around an n-fold axis), assuming that there are no correlations 

^ \ between the molecules. In the present paper molecular correlations are treated, but only in the sense 

. that several molecules could reorient simultaneously as in a cog-wheel mechanism. Moreover, we do 

' not examine the possibility that the relaxation times of these combined reorientations could depend 

CO , on details of the local environment created by the neighbouring molecules. Finally also an ergodicity 

■ condition has to be fulfilled. Admitting for all these assumptions we can show that the correlations 

. ^ ' do not affect the coherent quasielastic scattering pattern in the following sense: The functions of Q 
that intervene in the description of the intensities remain unaltered, while the functions of uj can 

I ' undergo a renormalization of the time scales. The latter changes cannot be detected as the time 

''7~| \ scales that would occur if the dynamics were independent are not available for comparison. In other 

, words: Coherent quasielastic neutron scattering is not able to betray the existence of correlations 

' of the restricted type that occur in our model. The assumptions that underly the model we present 
were made to allow a mathematically rigorous calculation of the scattering function. Other, perhaps 
more realistic cases may entail correlations of a type that is too difficult to be treated rigorously 

^ ' with our method of calculation. But our result presents an important non-trivial counterexample 

I , to show that the absence of a clue for the presence of correlations in the data is not a sufficient 

' O ' criterium to conclude that such correlations are indeed absent. 

■ 

O ' 61. 12. Ex, 87.64.-t, 66.30.Dn 

^ . I. INTRODUCTION 

O ; 

1^ i In a recent paper [1] we calculated the coherent quasielastic neutron-scattering signal [2-4] that results from n-fold 
' rotational jumps of N molecules periodically arranged on the sites of a d-dimensional solid-state lattice, by applying 
. a previously described method. [5] The jumping molecules themselves do not need to have n-fold symmetry. As a 
■ physical visualization of the present problem we could mention the array of nonadecane molecules hosted by an urea 
inclusion compound. [6] (Such molecules have CH3 end groups of three-fold symmetry, but the entire molecule has no 
' such symmetry). 

I In a first approach, we assumed that the molecules are jumping independently, leaving open the question: What 
H ■ could happen if the relaxational motions of neighbouring molecules cease to be uncorrelated? We anticipated that in 
I ' that case coherent quasielastic neutron scattering might be able to evidence possible correlations. In the present paper 
we address this possibility, understanding the word correlation in the sense that within connected sets of a given size 
and shape all molecules jump simultaneously in some phase or anti-phase relations, like in a cog-wheel mechanism. In 
Q , the paper we will use the word "clusters" to refer to these sets. Assuming further that the relaxation times for these 
combined jumps do not depend on the momentaneous orientations of the molecules inside (and outside) the cluster, 
we come to the astounding (and disappointing) conclusion that such correlations do not leave a fingerprint in the 
structure factors of the ensuing coherent dynamical signal. (See the abstract or below for a more precise formulation 
5—1 ■ of the restrictions). In other words: It is impossible to extract information on the presence and/or the nature of such 
correlations from coherent quasielastic neutron scattering signals alone. 

Let us point out where the surprise lies. Quite often, if we calculate the total structure factor for such dynamics 
within a single isolated cluster, we will be able to detect the presence of such correlations, as we do not find the 
same answer as for the jump dynamics of a single molecule. Even when the cluster is not isolated, highly correlated 
rotors still can be expected to manifest themselves very differently in reciprocal space than randomly independent 
rotors, e.g. in a linear chain anticorrelated configurations produced by the dynamics would manifest themselves at 
the Q- value corresponding to the zone boundary of the reciprocal lattice. Indeed, the correlated dynamics can give 
rise to static correlations between the orientations of the molecules, that should be observable. How is our result then 
possible? This question resumes one of the counter-intuitive aspects of our result. 
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The snag is that the clusters we consider can be centered on any point of the crystal lattice, such that two clusters 
situated on neighbouring points can overlap. If the cog-wheels of two such points have ground in succession then the 
rotations involved in the second move will partly undo the static correlations that existed between the molecules after 
the first move. In such kind of dynamics the correlations are no longer as all-out that we can call them "ferro" or 
"anti-ferro" (as in the example of anti-correlated configurations above). To complicate things further, they are not 
necessarily zero eighter. They can be of an intermediate strength. How the story unwinds further, really will depend 
on the specific dynamics defined within the cog-wheel mechanism. E.g. it is conceivable that in certain dynamical 
models the memory of the initial static correlations in a single cluster (or in the crystal), eventually are totally washed 
out with time, due to the game of overlapping clusters. We then have a different situation than the one evoked for a 
problem of an isolated cluster where the static correlations might persist indefinitely. Equally conceivable are other 
models, where such a total loss of memory does not occur. 

We will distinguish two mutually exclusive possibilities: (1) The loss of memory with time is total, such that starting 
from any initial configuration that one could obtain by independent dynamics we can stray to any other configuration 
that could be obtained by independent dynamics. We will call this the ergodic case. (2) On the contrary, in the 
non-ergodic case not all configurations can be reached. The configuration space that applied for the independent 
dynamics then splits up in several disjoint subsets. Only one of these is the set of possible states of the sample under 
study. All other states just do not exist in the sample. -"^ 

Neutron scattering calculations involve an average over all possible initial states. In the ergodic case, this averaging 
will run over the same set of initial states as in the case of independent dynamics. In the non-ergodic case, the set 
is smaller. We arc able to prove that this kind of averaging over all possible configurations of the molecules in the 
ergodic case is able to iron out the signature of a correlation that was present in the dynamics of an isolated cluster 
in reciprocal space, i.e. we recover the result for a single molecule. In the non-ergodic case, we average over a smaller 
set of possible states. This kind of weaker averaging does not permit us to assert that we would recover the result for 
a single molecule. That is the result of our paper. 

At this point it is perhaps instrumental to summarize the situation at the hand of an example. Let us take a 
linear chain whose 2N molecules can take only two orientations (which we can call pseudospins) . We suppose that 
the dynamics is given by some rule that two neighbour sites can flip simultaneously their pseudospins (by two-fold 
rotations). If this rule defines permanent dimers in the sense that only molecules 1-1-2, 3-1-4, 5-1-6, etc. can turn 
simultaneously, then we are in the presence of a strong correlation. The presence of such dimers is immediately 
obvious on an instaneous photograph of the system and will leave an unmistakeable fingerprint in the dynamics. That 
is the case we arc all familiar with. And it is on this case that we all built our intuition. 

We render the situation far more subtle by allowing the dimers to be dissociated by the game of overlapping 
clusters.^ Let us first consider a simultaneous flip of the molecules on sites 1 and 2 (starting from a completely 
aligned configuration). That introduces a correlation. Now assume that an analogous flip involving sites 2 and 3 
follows. Will this remove the "correlation" , as the "image" of the initial "dimer" has been destroyed? The answer is 
(here) that some less conspicuous correlation remains. Due to the construction of the model there will always be an 
even number of flipped pseudospins. How weak it ever might be, this is a (static) correlation that will never disappear 
from the model. In fact, due to the persistence of these static correlations, the system is not ergodic and does not fulfll 
the assumptions that are underlying our paper. But in this very specific case it happens that we can actually prove 
that the model does not yield the same neutron answer as a system with totally uncorrelated molecules. One can 
easily check this by using our method of calculation on a small toy model and exploiting the fact that the configuration 
space for the uncorrelated case splits into two disjoint subsets when we introduce the dynamical correlations. But 
if in another type of dynamical model the configuration space were not to split, then we will be able to assert with 



^Classifying the system according to this distinction was not motivated by the physics. This choice mainly reflects the limits 
of our abilities to perform certain calculations. The distinction wc make codes information about what happens to the system 
at infinite times: The correlations are eventually lost or otherwise. Nothing is stated about correlations that might persist 
for finite times. Despite this lack of information about the situation at finite times, in case (1) we will be able to infer an 
impossibility to discern indications for the presence of correlations within the dynamical signal. 

^There is then no longer equivalence between the concepts of a dimer and of a (dynamical) cluster as it is used by us in this 
paper. Wc use the terminology "dynamical cluster" to indicate the presence of a dynamical correlation, and to define sets 
of molecules that are supposed to be able to turn together, while the terminology "dimer" indicates the presence of a fixed 
orientational relationship between two molecules. Our use of the word "cluster" has also nothing to do with the idea of a 
domain wherein all molecules would have some special orientational relationships, e.g. to be aligned in a parallel fashion as in 
a ferromagnct. Our definition of the word cluster is a purely mathematical one. It is thus different from the usual intuitive 
physical idea. 
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certainty that the dynamical signal will not offer clues as to the presence of the dynamical correlations. 

II. MOTIVATION OF THE USE OF OUR METHOD 

A many reader might find the mathematical leap we will take to configuration space too cumbersome and too much 
of a fuss to his taste. Is all this really necessary? Part of the answer may reside in the following warning against a 
possible confusion. Consider the following (fallacious) argument. When we look at an instantaneous photograph of 
an ergodic system, the orientations of the molecules will in general look completely random. From this observation 
it takes only one step to argue that the "correlations are random". Would that not show that our result is trivial? 
[7] The problem is that the word "correlation" occurs with several different meanings in this context. The (spatial) 
static correlations between the orientations of the molecules that would be visible in an instantaneous photograph of 
the sample may well be random."^ However, the instant dynamical (or temporal) correlations are not random, as they 
tie up a well-defined set of molecules into a cog-wheel. 

To avoid confusion, it is thus important to realize that the word correlation can intervene with two entirely different 
meanings in the discussion and in the mind of the reader. We have static orientational correlations between the 
molecules of the whole crystal that we may detect in a snapshot of the system. We call such snapshots configurations. 
These static correlations arc the onc^s that matter in structural diffraction techniques. One will run headlong into 
trouble if one docs not properly sort out for oneself that moving to a dynamical context involves a shift of paradigm: 
We have to consider temporal correlations between pairs of such snapshots at all possible times. The cog-wheel 
mechanism defines dynamical correlations between molecules within a (small) cluster, in the sense that it tells you 
which molecules will turn simultaneously and how they will turn. The dynamical correlations are the ones that really 
count, not the static ones, although the former are of course responsible for the occurrence of the latter. In view of 
this remark, the loose phrasing "in the presence of correlations" appears thus as fraught with ambiguity.^ 

An important aspect of our method is that it develops a precise mathematical tool that permits to discuss all kinds 
of details that can arise in an argument about correlations in a very clear and rigorous way, without running into 
difficulties due to certain inadequacies of plain language. The examples of the correlations that persist on a linear 
chain despite the dissociation of the dimers and of the fallacy based on the ambiguous use of the phrasing "random 
correlation" as sketched above may serve to show how heedless it is to underestimate the difficulty of deciding on one's 
mere intuition if a correlation is present within a system or otherwise, and how discussing such matters in informal 
language can seed confusion. Two persons may easily use the same words and sentences to refer to mental pictures in 
their heads that are fundamentally different, without being aware of these differences. In a casual discussion it might 
e.g. take quite some time before it transpires that a person does not take properly into account the vital distinction 
between static and dynamical correlations we outlined above. The reverse of the medal is that mathematical rigor all 
too often leads to austere if not frustrated reading. 

III. POSSIBLE DOMAINS OF APPLICATION OF THIS WORK 

A. Scope of experimental probes 

We are presenting our work in terms of a coherent neutron scattering problem, but the formalism will also apply 
to other microscopic techniques, provided the scattering is really coherent and the technique allows to obtain a Q- 
dependence. Coherent scattering does not only occur with neutrons when the nuclei have an important coherent 
scattering cross-section. In contrast with neutron scattering, where the nuclei can have both a coherent and an 
incoherent scattering cross section. X-ray and light scattering processes jiggle the electrons and will uniquely lead to 
coherent scattering. Good candidates for applications are thus also quasielastic X-ray scattering and light scattering. 
In such techniques the Q-range tends to be small (Q » 0), but a lot of progress has been made in the last decade 



^In diffraction studies static correlations lead to diffuse scattering. [8] It is perhaps tempting to invert this implication in the 
sense that absence of static correlations implies absence of diffuse scattering, but this is just wrong in dynamics. 

In this language of static and dynamical correlations, our statement that wc do not consider environment-dependent dynamics, 
can be reformulated by saying that the jump times are not considered to depend on details in the static correlations between the 
orientations of the molecules inside or outside the dynamical cluster within a given configuration, but rather uniquely defined 
by the cog-wheel mechanism on this dynamical cluster. 
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in inelastic X-ray scattering techniques using synchrotron radiation, rendering available large Q-values combined 
with good energy resolution 1 mcV). [9] When there is no energy information, the energy-integrated quasielastic 
scattering will still yield the Q-dependence of the diffuse scattering.'' 

Extremely high resolutions can be obtained with the so-called speckle technique which allows one to study relaxation 
on a time scale of the order of 15 minutes! [10] Very good resolutions (neV) can also be obtained by nuclear resonance 
methods using synchrotron radiation [11] (e.g. exciting the 14.41 kcV Mossbauer level of ^^Fe), but in general this 
is an incoherent process (as in this process the radiation interacts with the nuclei rather than with the electrons). 
Under certain exeperimental conditions the process is nonetheless coherent, but then the quasielastic scattering has 
to be deteced against a very strong background of elastic scattering. 

B. Scope of experimental phenomena 

Coherent quasielastic neutron scattering studies are rare and this is certainly due to the difficulty of treating the 

problem theoretically. Incoherent quasielastic neutron scattering can often be modeled by a system of coupled rate 
equations, that one can solve analytically. But describing correlations, as probed by coherent quasielastic scattering, 
along such lines proves difficult. One solution to this problem is to drop the ideal of an analytical description 
all together and to make Monte Carlo [12] or molecular dynamics [13] simulations, but then fitting some variable 
parameters of the model to the experimental data will no longer be possible. 

There are two main categories of stochastic motion that lead to (non-magnetic) quasielastic scattering, viz. trans- 
lational diffusion and rotational relaxation.^ Our paper does not deal with translational diffusion.'' There exists a 
generalization of the Chudley-Eliott model that permits to calculate the coherent quasielastic neutron scattering signal 
for translational jump diffusion on a lattice. But this is an average description in that it is based on rate equations 
which rely on the assumption that the concentrations are rigorously homogeneous throughout the sample. [15] For the 
sake of completeness we refer the reader interested in theoretical aspects of translational diffusion to references [12] 
and [16]. The bulk of the experimental efforts have focused on the diffusion of D in NbDa; [17], Ag in silver halides 
[18] and in AgaSe [19], Rb in RbAg4l5 [20], Na in Na3P04 [21] Oxygen in UO2 [22] and F in fluorites [23]. 



^It may be noted that the formalism for inelastic X-ray scattering leads to expressions that link the scattering functions to 
the measured intensities in a way that looks identical to what we have for inelastic neutron scattering. However, important 
differences are very often passed under silence: The scattering functions that intervene in X-ray scattering are sums over all 
the electrons in the sample, while in neutron scattering they arc sums over all nuclei. But both are generally noted as S(CI,ljj), 
without any explicit reference to this important distinction. Thus some work has to be done to render the two different 
definitions of scattering function comparable. It is not a priori given that the electrons will just move as the nuclei in the 
photon field. In Raman scattering e.g. the polarizability of the atom intervenes and gives rise to a coupling factor that can 
be hard to evaluate and is different from the simple expression u • Q that occurs for neutron scattering. (Here u is the atomic 
displacement). The polarization of the electromagnetic radiation is important, etc... To our knowledge these points have never 
been worked out exhaustively for X-ray scattering. 

®In molecular liquids both occur simultaneously and they can be coupled. [14] 

'^Thcre is an unsuperable intrinsic difference between rotational and translational diffusion. The point is easily made by 
considering the translational jump diffusion of two particles on a linear chain (of n positions). If the two particles could 
both occupy a same lattice site simultaneously then the configuration space would be just a square lattice (whose points 
(l,l),(l,2),...,(i,j),...(n,n), correspond to configurations where the first particle is located at the site with coordinate i, and the 
second particle at the site with coordinate j). As real particles cannot be in the same site at the same time, the diagonal of 
the square lattice {i — j) represents forbidden configurations. Once these points have been removed from the square lattice, 
we are left with two totally disconnected mirror-symmetric triangles (the mirror-symmetry corresponds to identical particle 
exchange). Eighter of these triangles can be taken as the configuration space for the two particles. The generalization to TV 
particles leads to an A?^- dimensional simplex that is a (l/Ar!)-th slice of an AT-dimensional hypercube. The secular matrix that 
intervenes in the rate equations no longer has A'^-dimensional translational invariance in this truncated configuration space. 
The method of calculation we use in our paper is based on A'^-dimensional translational symmetry in configuration space and 
breaks thus down. (In models of correlated rotational jumps there are no such hyperplane cuts through configuration space 
corresponding to problems of excluded volume, but it is also a loss of global translational symmetry, splitting the configuration 
space in totally disconnected parts, that renders it impossible to perform the calculations for the non-ergodic case ). If it were 
possible to overcome such fundamental difficulties one could make a mathematically rigorous calculation of the total scattering 
function for a hard-sphere liquid, by taking a continuum limit of an appropriate lattice gas model. 
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What kinds of rotational relaxation is our paper then able to address? Let us immediately state that our model 
may not have a perfect counterpart in the real world.* It is a requirement of not breaking translational symmetry in 
configuration space that imposes the underlying assumptions of our paper. If these are not met, the calculation of 
the correlations is of a difficulty far beyond the methods we use. 

We noted at the beginning of the Introduction that the rotating molecules do not need to have the symmetry of 
the rotational jump. This is important. If the symmetry of the molecules is the same as the symmetry of the jump, 
then the sample will look the same after the jump as before the jump. This means that the Fourier transform of the 
molecules can be factorized out in the calculations and we do not need to go through the more difficult approach of 
our paper. ^. 

Combining these two limitations of feasibility and technical interest, the geometry of the nonadecane molecules is 
intrumental for visualizing the problem that we can treat, and for describing the (restricted) kind of correlations that 
we have in mind. The physical correlations between the nonadecane molecules that arc really experimentally observed 
in the urea inclusion compound are of a different nature [27], which wc arc unable to treat (although this was not 
yet clear at the time we started this work). But this kind of correlations could exist in another system. The major 
lesson from our exercise must be that we give a non-trivial counter-example that indicates how coherent quasielastic 
neutron scattering may not always yield the information one expects on the basis of commonly accepted notions. 

In our method the configuration space for the single-molecule dynamics is one-dimensional. It is however possible 
to generalize the formalism such as to recover also higher-dimensional single- molecule configuration spaces. ^° This 
broadens the scope of possible applications further to three-dimensional rotations of molecules whose symmetry is not 
the same as that of the rotation symmetry. An example of such cases would be the rotational jumps of molecules of 
tetrahedral symmetry containing different isotopes (such that the symmetry is broken) as e.g. CH3D, CD3H, CH2D2, 
NH3D, C^'^Brs^iBr, C^^Bra^^Bra [24], CBraCla [28], CBrCls [28], etc... However in the physical world these molecules 
do not necessarily give rise to solids where our type of problems arc encountered. In partially deutcratcd methane the 
physics is just different: Bulk CII3D and CD3H are either close to rotational diffusion, or there is quantum behaviour. 
In partially deutcratcd ammonia, the existence of exchange of D and H between the molecules complicates the picture. 

Also octahedral molecules could be considered. Finally, an important application of our result could be Ceo- 
The (coherent) quasielastic neutron scattering signal of the rotational dynamics of this molecular crystal in its room 
temperature phase has been successfully described in terms of a model of completely independent rotors by Neumann 
et al. [29]. Although our results have been derived for rotational jumps rather than for continuous rotations as occur 
in Ceo, they suggest by analogy that it cannot be claimed with certainty on the basis of the sole neutron data of 
Neumann et al. that the molecules would be totally independent. 

IV. METHOD 

The general method we use has been already described previously [1,5]. The idea is to formulate the problem in 

configuration space. The instantaneous configuration of the whole system is an abstract particle that diffuses on a 
network in this configuration space. The vertices of the network are the possible configurations of the system. (How 
to define them is illustrated by our example of the translational diffusion of two particles in footnote 7). When one 
single particle of the system makes a jump with a relaxation time r, we say that the system "jumps" between two 
configurations with the relaxation time r. Two configurations that are linked by such a move of one particle are 
connected by a line of the network carrying a label r. The dynamics of our system is this way mapped isomorphically 



Eventually the assumptions that underly our model are dictated by the limits of our mathematical skills rather than the 
physics we attempted to address. As we will see, our approach of the problem leads to a description of the dynamics in terms 
of huge jump matrices. Analytical diagonalization of matrices of rank larger than 4 is in general not possible: As Galois has 
shown, beyond degree 4 there are no general solutions in terms of radicals for polynomial equations. A notable simplification 
can occur however, if we cam make use of symmetry. It is mainly translational symmetry in configuration space that achieves 
the analytical diagonalization of the jump matrix in the present paper. 

^Unless the dynamics is more complicated than a simple rotational jump model, due to supplementary problems such as 
translation-rotation coupling [24-26] etc... 

'^'^We will not give this derivation in order not to burden the paper. The proof runs as for the one-dimensional case mutatis 
mutandis. 

^^A good example would bo SFs, if it were not that the natural element fluorine is mono-isotopic (^^F), such that the 
simplification mentioned above can be applied. [25] 
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onto the problem of the diffusion of a single abstract particle (the system) on a network in configuration space. The 
only difference is that the embedding configuration space is of considerably higher dimension than physical space. 

The physical problem can thus be described in terms of a set of coupled rate equations -^P = MP, where the 
column vector P{t) contains the probabilities Pj{t) that the system is in configuration j at time r is a relaxation time, 
and M is the so-called jump matrix. The formal solution is worked out in reference [1] and does not contain further 
suprising steps or new ideas: By diagonalizing M = SAS^^, we can write the solution of this set of coupled linear 
differential equations as P{t) = S exp (Ai/r) P(0). By plugging in the various initial conditions into P(0), we 
obtain then the full set of probabilities Pj,k{t) that the system is in configuration j at time t if it was in configuration 
k at time 0. These are the quantities that arc needed in the Van Hove correlation functions. A proper choice of 
initial conditions will also take into account the thermal occupation factors, i.e. entail the required thermal averaging. 
The total scattering function is obtained from the spatio-temporal Fourier transform of the thermal average of the 
correlation functions. In its final expression the total scattering function reads: 5(Q, uj) = -hq- F S JF( e'^*/'^ ) F^, 
where F is the row matrix that contains the spatial Fourier transforms Fc of the configurations c obtained by putting 
a Dirac measure of weight bx at the position of each atom of type x {bx is its coherent scattering length); stands 
for the temporal Fourier transform. The normalization factor -\- is here written for the case of uncorrelated n-fold 
rotational jumps of N molecules. 



V. FORMULATION OF THE JUMP MATRIX 



When there are no correlations between the molecules, the configuration space for the n-fold rotational jumps of 
N molecules on a lattice C C R'^ {d e {1,2,3}), e.g. C ={[!,£] CiN)'^, N = i'^, will be a hypercubic lattice [30] 
H = {"L/n)^ with cyclic boundary conditions in R^. (Each particle adds one dimension to configuration space: E.g. 
if there were just two molecules (A'' = 2) that make threefold rotational jumps (n = 3), the configuration space would 
be the square lattice of the = 3^ = 9 points (1,1), (l,2),...(i,j),... (3,3). The configuration corresponds then 
to the situation that particle 1 has orientation i and particle 2 orientation j. The boundary conditions are cyclic as 
each particle can go through the successive orientations 1 — > 2 ^ 3 ^ 1). For a three-dimensional sample (d = 3), 
the position vectors of the molecules will be {jx,jy,jz ) S C The jump matrix M is then defined by: 

Mc;d = -2A4;d + Yl ( ^•=;d-Kej + (5c;d-ej )• (1) 

iec 

The reader should not feel intmidated by the few very concise notations of this type that occur in the paper. They 
are mainly givc;n for the sake of completeness. The main point he should capture is at which points in the jump 
matrix there are non-zero entries and how they look like. Eq. (1) is built up as follows. There is a unit vector ej 
in hypercubic space associated with each molecule (with position vector j) that can turn: As we illustrated with 
the example of two particles above, each molecule adds a dimension to configuration space. When we start from a 
configuration c and the molecule j makes a rotational jump, the system goes to the configuration d = c -|- ej. Hence, 
in the absence of correlations, a configuration c e W has 2N neighbours d = ciej, where j = {jx,jy,jz ) G £: At the 
line corresponding to c in the jump matrix there are thus non-zero entries at all columnar positions corresponding to 
the configurations d = c -|- ej and d = c — ej, a fact which is expressed through the presence of the Kronecker delta's 
^c;d-i-ej and (5c;d-ej • The diagonal terms follow from these terms as in any other jump matrix for a diffusion problem. 



^^Generally spealcing, treating the system in terms of configurations has several advantages. First of all, the total scattering 
function is invariant under the operation of identical-particle exchange, and this fact is automatically accounted for in our 
method, by considering as strictly identical all the configurations that can be obtained one from another by swapping identical 
particles. Consequently, identical particles are not tagged in our method: They remain undistinguishable. We do not loose 
time or efforts in enumerating tagged configurations. This reduces the number of configurations to be considered drastically: 
For each type of particles j, there can be a reduction by a factor nj\, where rij stands for the number of such identical particles 
of type j (provided that they can indeed physically swapped by the dynamics). This is e.g. the reason why in the example of 
the translational diffusion of two particles in footnote 7, the configuration space can be taken as only one of the two triangles 
defined by the diagonal. (In the present case of rotational jumps, this simplification does not occur as the dynamics do not 
entail exchanges of molecules). The second advantage is that all static correlations that might occur inside a single configuration 
are automatically taken into account. In the calculation this transpires through the presence of the Fourier transform Tc of the 
configuration c. For this reason we further will never even mention these static correlations in the paper. The worth correlation 
only occurs in a dynamical sense. 
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The unit vectors ej are the generators of the hypcrcubic lattice and each configuration c accessible can be 
written as c = cj ej, with cj G Z/n. The hypercubic lattice Ti represents all possible configurations that can be 
obtained by n-fold rotations, and therefore the configuration space when there are correlations will be a sublattice 
S CH oiU}^ 

It is logical that we assume that due to the range of the intermolecular interactions, the set or cluster of jumping 
molecules extends to a neighbour shell of a certain order. Due to the translational invariance on the physical lattice 
C, this cluster should be allowed to occur at every lattice site j e £. A configuration c will now be connected to other 
configurations c ± vj, j e £ on the sublattice S, than c ± ej, j G £ as for the independent dynamics, and it is the set 
of V relative position vectors vj that has to be defined. E.g. if on a two-dimensional square lattice C ~ ( fl N)^ 
with cyclic boundary conditions, the jump of a molecule at {jx,jy) over 27r/n is always accompanied by opposite 
jumps over — 27r/n of its four first neighbours, then v = i'^ and 

= ^ij.,3y) -^ij.+l,3y) -^(j.-hjy) -^U.Jy + l) - ^ U . ,3y -1) i J X , jy ) ^ A (2) 

where the minus signs translate the fact that the rotations are opposite. We can understand this by decomposing the 
combined move mentally into a succession of single-molecule jumps: We first turn the molecule {jx,jy), which changes 
the configuration from c to c + e ^^ j^^y Then we turn the molecule {jx + ^, jy), which changes the configuration further 
from c + eQ-^ j^) to c + e{j.^.,jy) — ^{j^+i.jy)^ etc... upto the final configuration c + '^(^j^j^). The order in which we 
take these individual moves is not important, as in the combined move of the five molecules the system does not visit 
the intermediate configurations and goes immediately from the initial (c) to the final configuration (c -|- v j ) ) . On 
our network there will be a line that connects these two configurations and which will be labeled by r. o/ course 
the relaxation tim,e t we use now has no longer any relationship whatsoever with the relaxation times we used in the 
problem of independent dynamics. It is the relaxation time for a cluster, not for a molecule, and the probabilities we 
are dealing with here can in principle not be expressed as a product of probabilities as would occur if the rotors were 
independent. ^''^^ The v = vectors Vj,j = (1,1), (1,2), . . . {£,£) are the generators of the sublattice <S. Using the 
hypercubic norm, we see that || '^(j^jy) II = 5. 

As a second example, we could also imagine that each molecule at (jx , jy) has (for symmetry reasons) four equivalent 
alternatives to jump simultaneously with a single first neighbour. Then z/ = 4^^ and y{jx,jy) G 

^i3.,3y),3 =^U.,3y) - ^OxJa-M)' 
^Ux,jy)A — ^UxJy) ~ ^Uxdy-i)- 

Here || 'V(j^j^)^r 11= 2 ( r G {1,2,3,4} ). We see thus that in the general situation with correlations the jump matrix 
M will be defined by 



We would like to stress that there is no underlying assumption related to representing the various configurations as points 
on a hypercubic lattice, for the mere convenience of enumerating them more easily. This also true when the rotors are not 
independent. More specifically, it does not imply any factorization of probabilities as we would have in the case of independent 
rotors. The probabilities will be given by the topological connectivity of the network or graph that links these points. These 
lines could be very different from the edges of the hypcrcubcs of the hypcrcubic lattice that occur in the case of independent 
dynamics. The probability for a jump between two configurations is taken into account by connecting the points that represent 
them by a line labeled with its relaxation time. Only if these probabilities themselves were factorized would we have an 
underlying independence. 

^*It must be obvious that the occurrence of a cluster is probabilistic: The model of Eq. (2) is perhaps best seen as part of 
a larger model wherein e.g. a molecule has some probability to turn alone, some probability to turn as given by Eq. (2), and 
some probabilities to turn even within other types of clusters (see e.g. Eq. (13) below). One should thus not use the definition 
of Eq. (2) to "prove" that the whole lattice should turn simultaneously: when {jx,jy) turns, {jx + 'i-,jy) turns; then since 
{jx + 1, jy) turns, {jx + 2,jy) must also turn, etc... That would be an absurd way of misinterpreting of our aims. 

'^^We must emphasize that further on our approach will only work since we assume translational invariance in configuration 
space Ti: the types of possible rotational jumps of a molecule (or a cluster of molecules) and their respective probabilities do 
not depend on the specific orientation of that molecule or any other actual detail of the prevailing configuration of the system. 
This would be quite reasonable for the case of alkane molecules evoked, but nevertheless this remark implies that we have 
not exhausted all possible meanings of the concept correlation. (We arc not dealing in the present paper with the additional 
complication that the dynamics could also depend on the actual local environment of the molecules. Defining systems of the 
latter type in a self-consistent way is already difficult in its own right. We briefly touch upon it in the Appendix.). 
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p 

Mc;d = -2Np5c-A + ( '^<=;d+Vj,. + -^Cid-Vj,, ), 



(4) 



where p = v j l'^ is the number of difFcrcnt types r of clusters per lattice site that can turn (i.e. p = 4 in the example 
of Eq. (3)). The definition of Eq. (4) is analogous to the one in Eq. (1) with each vector ej replaced by p vectors Vj,r- 
We can more generally assume that the various possibilities labeled by r lead to different time constants. Drawing in 
the time constants into M the jump matrix becomes: 

P 1 

Mc;d = XI ~ ( -2^^c;d + X ( '^c;d+vj,. + <5c;d-vj,. ) )• (5) 
r-=l je£ 

This notation can be further generalized by noticing that r corresponds to a difference (or translation) vector t € £, 
l^te£- Wc must then admit that for many values of t, we will have = 0, or replace t G £ by t e ^ where Q 
is a physical cluster: 

Mc;d = X - ( -2^^c;d + X ( ^c;d+vj.t + -^cjd-vj.t ) )• (6) 

tee '^^ je£ 

We have just given some examples of possible correlations. In the present paper we will not further insist on writing 
down a general formalism that would cover all possible cases. First of all, there are too many possibilities. Secondly, 
we think that a general abstract formalism in configuration space [31] (which typically has a dimension w 10^^) with 
notations of the rather elaborate type that occur in Eq. (4-6), might easily conceal the rather simple ideas behind 
our method. We will therefore rather proceed by illustrating the method for examples of the type given, in order to 
familiarize us with its spirit. 

VI. THE PROBLEM OF ERGODICITY 

Since the hypercubic norms of the vectors vj are by definition larger than 1 in the presence of correlations, one 
would be inclined to think that <S will always be a strict subset of H, which we expressed by saying that our dynamical 
problem is no longer "ergodic" . The real situation is more subtle. Taking the possibility of non-ergodicity seriously 
is a difficult task that we arc unable to treat rigorously in general. In the Appendix we develop a physical argument 
to put forward the idea that we can relax the condition of ergodicity to a weaker criterium of local (as opposed to 
global) ergodicity. For the main stream of the paper we will henceforth assume ergodicity. But despite the fact that 
we devote here only a few lines to this assumption, we must stress that ergodicity is a vital issue for the validity of 
the calculations in the rest of the paper. 

VII. SOLUTION OF THE JUMP MODEL IN THE ERGODIC CASE 

When the system is ergodic, we can immediately use the eigenvectors we already established for Ti in reference [1] , 
and we will be able to derive without effort the eigenvalues from the form of Vj^^, j S C Indeed, when the system is 
ergodic then the jump matrix in configuration space has translational symmetry along each of the directions defined by 
the N unit vectors ej that span H. The eigenvectors are therefore just A^-dimensional Bloch waves, i.e. x 1 column 
matrices V^'') defined by [ VC') ]c = cxp [i ^ (k - ko) • (c - ko) ], Vc G U. Here ko stands for (1, 1, 1, ... 1, 1) e 
and is introduced to take into account the fact that the expression for the eigenvectors features the quantities fc^ and 

always under the form of linear combinations [k^ — 1), and (cj — 1). The iV-dimensional Bloch wave is obtained 
as a Kronecker product of N one-dimensional Bloch waves, just as a three-dimensional phonon can also be written 
as a Kronecker product of three one-dimensional Bloch-waves. (Each one-dimensional Bloch wave is a n x 1 column 
matrix). Combining these eigenvectors with the definition of the x jump matrix M in Eq. (4) yields the 
corresponding eigenvalues: 

p 

Ak = -2p7V + 2 X X^°^[^-i'--('*-'*o)J- (7) 

This is exactly analogous to the way the eigenvalues of a phonon problem are obtained by operating the dynamical 
matrix on the Bloch wave eigenvectors. We just calculate MV^"*^ using the definitions of M and V^"*^ and check that 
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the result can be rewritten as 

A(k) yi^) for sonic value A'^*'^. For the phonon case in text books, this calculation is 
usually only written for a general line of the column matrix that represents the eigen vector, with the possible effect 
that the argument might not be recognized as perfectly analogous to the one we are dealing with here. E.g. in the 
case of Eq. (2) we have A'' = p = 1, Vj € £ : ko • vj = —3, such that Eq. (7) is more explicitly seen to yield 

Ak = -2^2 +2cos[^(fc(i,i)-fc(i,2)-fc(i,^)-A;(2,i)-A;(^,i)+3)] 

+2 cos [^( -A; (1,1) +fc(i,2) -fc(i,3) -fc(2,2) -fc(€,2) +3)] 

; (8) 
+2cos[^(fcQ-^j^) -fcQ-^+ijj -fcQ-^_ij^) -A;(j-^j^+i) -fc(j^j^_i) +3)] 

+2cos[^ - kf^i^e) - k(e-i,e) - k(e,i) - k^e,e-i) + + 3)], 

where k — ( fc(i i) , fc(i,2) ; • • ■ ■ • • ); and k{j^,jy) € Z/n. There are thus n^^ k-vectors in the case of 

Eq. (2). In fact, a general configuration c has £^ vector components: c = (c(i,i), C(i,2)) ' ' " C(j^,jj,), • • • C(^,^)). The 
coordinates cq^j^^ can take n values representing the n possible orientations of the molecule at {jx,jy) G ^- There 

are thus in total n^^ configurations, and for each of them there is a probability. The jump matrix M works on the 

space of these configurations. An eigenvector V'^'') is a vector with entries [V'^'^^Jc, where c runs through 
all possible configurations. We can thus depicture V''') as a function (the Bloch wave): V^^^) : c [V('')]c. Let 
us operate M on an eigenvector. The first non-zero off-diagonal entry in M corresponds to the possibility that c 
undergoes a transformation to c + v (i i) . That will transform the function from exp 221 (k — ko) • (c — ko) ], Vc e W 

toexp[z^ (k-ko) • (c + V(i_i) -ko)],Vc € H, i.e. multiply it by exp[i^(fc(i,i) -A;(i,2) -^(i,^) -fc(2,i) -fc(€,i) + 3)]. 
Continuing this way with the possible transformations towards c + V(i_2),-") c -|- y{e,e), including also the possible 
transformations towards c — V(i_i), c — V(i,2)v) c — V(^_^), and finally adding the diagonal term we obtain the result 
announced in Eq. (8). 

The whole part of the calculation involving the determination of the structure factors remains the same as for the 
case without correlations. In fact this calculation (based on the evaluation of G = F S) is entirely defineded by the 
values of the eigenvectors (S) and the Fourier transforms of the configurations (F) and these are not changed in the 
ergodic case. This means that Eqs. (35) and (37) of reference [1] remain valid. Using the same notations as in that 
paper, the Lorentzian L(hX-k_/T, cu) will be associated with | Gk | ^, where 

Here ^c(j^ is the Fourier transform of the molecule {j^ , jy) (with orientation cy^j^) ) if it were placed at the origin. 
(The real position vector of this molecule is T^{j^,jy))- Consequently, the structure factor of the elastic term remains 
the same. The structure factors of the quasiclastic terms remain also the same. Just the widths of the Lorcntzians 
have to be re-examined. Due to the presence of the Kronecker symbols in Eq. (9) only Lorentzians for k-values 
which have one single component k^^j^^ ^ 1 are contributing. In our example on the square lattice, these are the 
Lorentzians with 

27r 

Ak = -2^2 + 2 {i^ - 5) + 10 cos [ — ( A; 1 )] , (10) 

since in Ak there arc — 5 terms with a vanishing argument in the cosines and 5 terms where the argument of the 
cosines collapse to ^(^(j^j^) — !)• We recover thus miraculously the same Lorentzians with a width parameter 
—2 + 2 cos [^{k(j^jy) — 1 ) ] as in the case without correlations, except for the detail that the relaxation time r that 
occurs in the rate equations is being replaced in the final result by a five times faster relaxation time t/5, which does 
not happen in the case without correlations.^^ There is only a possible renormalization of the relaxation time to be 



But since we are dealing with clusters of 5 molecules, we might originally have had the intuition to put a value 5 times 
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used in the description of the jump dynamics. Conclusion: It is absolutely impossible to appreciate from a signal due 
to coherent neutron scattering alone if there are correlations or otherwise! One can only hope that a comparison with 
incoherent data might still yield some clues. This result will be general if the cardinal number of the set of possible 
configurations H is not altered by the presence of correlations (case (1) in the terminology of the Introduction). In 
fact, the factor 5 in our example comes from the circumstance that our clusters contain five molecules, such that 
^Ox-jj,) ^'^11 contain five terms. Due to this Ak will contain five terms that contain ( ^ 1 ), each time combined 

with other terms of the general type ±( ^(jv ) — 1 ), that vanish since ) = 1. All other contributions to Ak are 

devoid of terms in k(^j^j^y 

We have not addressed the possibility illustrated by the example of Eq. (3), where we have more than one type of 
cluster. In principle this case is not different in nature from the one embodied by Eq. (2), since there we also have 
already u > N, such that Eq. (2) represents two equations rather than one. Provided we have ergodicity we can 
derive for this model from Eq. (4) with A/' = p = 4, (Vj e £)(Vr e {1,2,3,4}) (ko-vj,^ = 0): 



Ak = -8^2 -|.2cos[ 




-fc(l,2))] 


+ 2 cos [ 


:f(^(M) 




+2 cos [ 


V(*(M) 


- fc(2,l) )] 


+ 2 COS [ 


:v(fc(M) 


" ^(£4) ) 


+2 cos [ 




-^(1,1))] 


+ 2 COS [ 


:f(fc(i,2) 


-^(1,3)) 


+2 cos [ 




-k(2,2))] 


+ 2 COS [ 


:f(fc(i,2) 


- fc(«,2) ) 



(11) 

+2 COS [ f ( A; (,-^ ,,-^) - k ) ] + 2 cos [ f ( A: (,-^ - k (,-^ ) ] 
+2 COS [ ^ ( ,,-^) - k ,,-^+1) ) ] + 2 cos [ ^ ( ,,-^) - k ) ] 

+2 cos [ ^ ( A; (^,^) - fc (1 ,^) ) ] + 2 cos [ ^ ( A; (^,^) - fc 1 ) ] 

+2cos[^(A;(^,^) - A;(^,i))] +2cos[^(fc(^,^) - A;(^,^_i))], 

Expressing the selection rule from Eq. (9) that only Lorentzians are contributing which correspond to a k-value that 
contains a single component A(jj,,j^) 7^ 1, we just retain the Lorentzians with 

27r 

Ak = -8^2 + 2{4:e-^) + W cos [—{k^j^j^^-l)], (12) 

such that the widths are now again —2 + 2 cos [ v ( ^ (ix jj,) ~ ^ ) ] ' this time with a prefactor 8, since there are in total 
16 different moves that involve a jump of molecule {jx,jy )■ 

We may note that we have already covered quite realistically looking possibilities. The task of taking into account 
correlations of this type on a whole lattice definitely seemed daunting before we embarked on our method. The most 
realistic case would probably involve a distribution of cluster sizes, whereby each cluster has its own characteristic 
relaxation time. E.g. on the square lattice we could instead of Eq. (2), have relaxation times and clusters y{jx,jy ) G 
£: 

"^1 • ^Ux,3y),i = -e(j,+ijj -eQ-^_ijj -e^^^^+i) -e^^j^^i), 

'^2 : '^Ux,3y),2 = ^U.dy) -^U. + lJy) -^U.-lJy) -^{j.dy + l) -^{j.dy-l) + « (j, - 1, +1) + 6 (j^ + 1 _ 1) 

+ ^{jx+2,jy)+e{j^-i,jy-i) +eo,-2,j„) +e(j^j^+2) +e(j-^j^_2) +e(j^+ij^+i). 



slower than in the uncorrelated case into the rate equations. If we do not allow for this, then we would end up with a speed-up 
of the quasielastic signal in the presence of reorienting clusters, which is quite unphysical, as it should in principle be more 
difficult to move a set of molecules rather than a single molecule. This use of a factor 5 based on a handwaving argument can 
only be considered as a first approximation. We have no means to elucidate this relationship any further, since the jump time 
a single molecule would have had if its jumps had been independent is just not available to the experimentalist. There could 
be a change in the jump time parameter between the cases of correlated and uncorrelated jumps, but this cannot be learned 
from an inspection of coherent scattering data alone, since there is no change in the structure factors. 
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(13) 



whore in the additional cluster we have included now all members of the second-neighbour shell, assuming that their 
rotional jumps are in phase with those of the central molecule, i.e. in opposite phase with all the members of the 
first-neighbour shell. This leads to Lorentzians with a width parameter — A^^'' -|- — Aj^^' , where A|^^'' is given by Eq. 
(8), while A^f^ = -2^ + 2 Eq-, e£ ^oa [ ^ {k^J^^J^)-k(J^+^^J^)-k(J^_^^J^)-k(^^^J^+t)-k^^J^^J^_r)+k(J^_^^^^+l) + 

+ ku^+2,jy) + k(j,-uj^-i) + k{j^-2,jy) + fc(i.,j„+2) + fc(j.,j^-2) + A;(j^+ij^+i) - 5)]. We can appreciate 
that in general the width parameter for the Lorentzian Lk will be 

^l(_2iV + 2^cos[^Pj,,]), (14) 

since each cluster, i.e. each term vj^^. introduces a term — 2A'' + 2 Ej g c ( ^''^ ^j,^/^ )> where the polynomial Vj^r 
in the fc's is obtained by replacing each symbol e by a symbol k in the definition of vj^^ and adding — ko • vj After 
combining with (5ifc (11) ^ife(i,2) ••• '^ifeo^j^-i) ^ifco^.j^+i) ••• "^ifew,*-!) <5i fc in the example of Eq. (13) we recover 
— 4( + sin^ [^{k — 1 ) ] provided S = H. (Here we have used the fact that the set of values k(^j^j^^ G Z/n is 
the same for all molecules {jx,iy), such that identical values coming from different values of (jxTjy) can be regrouped 
in a term where the A;- value is simply noted as k). The more different types of clusters one includes, the more likely 
the configuration space <S will cover the whole of H. In fact, if also single-molecule jumps are allowed we forcedly 
end up with S = H. Also for the final expression of the Lorentzian widths after the application of the selection 
rules it is not difficult to derive the general result; —4 [ Er=i ( "^r/^r ) ] sin^ [ ^( /c — 1 ) ], where rrir is the number of 
molecules involved in cluster r. The use of the quantity rrir exploits the fact that in general the correlations between 
the molecules will be symmetrical. The most striking point is that all this follows almost effortless, even for very 
complex situations, as the translational invariance on H establishes a kind of dictionary vj^^ 'P},r which allows 
immediately to write down the eigenvalues, and that these eigenvalues are further tremendously simplified by the 
selection rule for the k-values. 



VIII. CONCLUSION 

We conclude: Apart from a possible renormalization of the jump time, the signal in presence of correlations of the 
type we considered is the same as in total absence of correlations. In other words: Coherent quasielastic neutron 
scattering is unable to reveal possible correlations of the type we considered between the rotational jumps of the 
molecules. This is a clear, but non-trivial result. We may call it paradoxical that coherent scattering cannot reveal 
the type of correlations we considered, while incoherent scattering can. The finding that the widths are renormalized 
calls for caution in the interpretation of such widths in the results of an experiment. One would be inclined to interpret 
them as simple quantities corresponding to the relaxation of a single molecule, while their real meaning could be quite 
different. That is another important result of our work. 

Acknowledgements. The author would like to thank M. Duneau for fruitful discussions. 

APPENDIX: THE PROBLEM OF ERGODICITY 

A. Why the problem is not simple 

As the hypercubic norms of the vectors vj are by definition larger than 1 in the presence of correlations, one might 
believe that we never have ergodicity. They span a larger fundamental cell in hyperspace than the unit cell that 
occurs with independent dynamics, and therefore less of such cells should be contained in H. But since we have cyclic 
boundary conditions in each dimension of the hypercubic lattice H this is not necessarily true. A good example to 
illustrate this point is the transformation vi = e^; + e^, V2 = e^, — on a square lattice, which splits into two 
sublattices, one with integer and one with half- integer coordinates in the not normalized basis vi, V2. But if we apply 
this transformation to the restriction (Z/n)^, allowing for cyclic boundary conditions, then for n odd the vectors ei 
and 62 will generate the whole lattice. We see thus that the problem of ergodicity we ran into is not as trivial to cope 
with as one might have thought at the start. It is however very important to know if we can claim that the problem 
is ergodic, since it determines if the calculation presented in this paper applies or otherwise. 



11 



Let us consider in the non-ergodic case a configuration Ci that occurs in 71, but not in the sublattice S. On this 
configuration we can apply also the set of generators Vp,p = 1, 2, ... i^, to obtain the coset Ci + S. The term coset is 
here appropriate since the lattices Ti and <S can be considered as visualizations of a translation group. There might 
be a further configuration C2 that does not occur within <S U (ci + <S). This way we construct s = ^Ti/^S disjoint 
cosets, based on initial configurations ci, C2, . . . c^. Our lattice H decomposes then as a Bravais lattice S convoluted 
with a unit cell containing s "atoms" Ci, C2, . . . Cg. 

The correlated dynamics will then map onto the diffusion of an abstract particle (the configuration of the system) on 
the connected graph whose vertices are the points of S. From any point of this graph it will be possible to wander to 
any other point of this graph. A point from another coset cannot be reached. Another coset represents thus a similar 
sample with similar dynamics, but the states of the sample represented by this coset cannot be reached starting from 
the states of the sample represented by S. The s cosets represent thus s different samples with a similar structure 
and similar dynamics, that would be a kind of "eniantomorphs" . The different samples arc characterized by the fact 
that they cannot have the same initial state. It is a priori not true that working with one of these s enantiomorphs 
should still yield the same Q-dependence as for the independent dynamics. This is shown by the counter example 
with the Ising pseudospins given in the Introduction. 

The main problem is that we would have to do the thermal averaging over S rather than over Ti. We would therefore 
need to determine how the lattice S looks like and to find a set of basic vectors that generates it. In general, this will 
prove just too difficult a task, due to the huge dimensions of the configuration space. Furthermore, the averaging over 
S will in general be too weak. Many sums in the averaging over H reduce to Kronecker deltas (see reference [1]). In 
the average over S, these sums will contain less terms, and therefore no longer necessarily yield results that are either 
or 1. 



B. Obvious mathematical criterium of non-ergodicity 

One can think of one simple way to find out that the sample is not ergodic, viz. when the dynamics can be described 
by a single equation, as e.g. in the example of Eq. (2) and the dimension of S is lower than the dimension of Ti.. 
There is a linear transformation with matrix T: 

[V(l,l), V(i_2), • ■ • V(2,l), . . V = T [e(i_i), 6(1^2), . . -6(2,1), • • • ^ ]^ ■ (15) 

If det T = 0, it is clear that the dimension of S is smaller than the dimension of H. It is instrumental to note that 
— (T + 3 H) is upto the factor 1/t equal to the jump matrix A for the translational diffusion on the square lattice 
£ = ([1,£] n N)^, which we know to diagonalize due to its translational symmetry when we assume cyclic boundary 
conditions. For other types of correlations than given by Eq. (2) an analogous argument will hold. The eigenvalues of A 
are -4sin2 [ f ( fc^-l ) ]-4sin^ [ f ( fcj^-l ) ], hence the eigenvalues of T are -3+4sin^ [ f ( ^^,-1 ) ]+4sin^ [ f (fcy-1)], 
such that: 

i £ 

detT=n n {-3 + 4sin2[f (fcx-l)]+4sin2[^(/e,-l)]}. (16) 

If one of the factors that occurs in this product is zero, we will know that the problem is non-ergodic. E.g. for 
£ = 12 & fca; = 3 & fcy = 4, we have —3 + 4 sin^ [ j ( fca: — 1 ) ] + 4 sin^ [ f ( ~ 1 ) ] =0, such that in all cases where 12 is 
a factor of £, the determinant will be zero, leading to a reduction of the dimensionality of the problem. We can infer 
from this that our dynamics very often will not be ergodic when we adopt cyclic boundary conditions on C We will 
therefore sketch a non-rigorous argument that suggests that for large samples non-ergodicity is not important. 

C. A loophole of escape 

Let us imagine that the neutron beam only illuminates a small part (£) of the whole sample (£'). We have then 
C C C. We assume that the dynamics on C are not ergodic. Let us call H' the configuration space that can be 
constructed from C and Tii the configuration space that could be constructed by independent dynamics from £, if it 
were an isolated piece of sample. Let us drop from the description of the configurations in 7i' everything that cannot 
be seen by the neutron beam, i.e. restrict their description to C. Very often this incomplete, truncated description 
will generate all possible configurations of C, i.e. the whole of Tii. In other words, the system is then not ergodic but 
the neutron beam cannot see it: £ has a (fake) apparent ergodicity. This boils down to a weakening of the ergodicity 
from a global to a local criterium. On any sizeable patch of the sample all possible configurations will then occur. 
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We will be in such a case of local ergodicity if we can build from a linear combination with integer coefRcients of 
vectors Vk, a configuration cj G 7i', whose restriction to L is equivalent to ej, for some j G put differently, if the 
projection of this configuration onto Jii is ej. By translational invariance, we will then be able to do this Vj G £, and 
thus we will also be able to construct by linear combination any configuration c gHi, i.e. we can build a set whose 
projection is the whole of H. 

In the case of Eq. (2) it looks quite plausible that this would be possible: We start with vj. By using Vk 
contributions from the second- neighbour shell we can zero the coefficients of the ek corresponding to the first- 
neighbour shell, and so we can proceed further outwards shell after shell, upto arbitrary distances from j. For the 
example of the dynamics given by Eq. (2) we have verified in detail that this works out. There are no conclusions to 
be drawn from this single case for any other posibble model. One will have to check case per case if one can fall back 
onto such an argument of local ergodicity. 



IX. APPENDIX 2. DEFINING ENVIRONMENT-DEPENDENT DYNAMICS 



We want to show here how one can define environment-dependent dynamics solf-consistently. This will allow us to 
explain why we are not able to treat such dynamics by our method. Let c € H be a configuration, c J2jec ^J- 
we want to jump from c to c ± ej we must look at the environment of the molecule at j, i.e. we must look at the local 
configuration built by all the molecules that belong to a cluster j + G centered around j, where Q is a set of relative 
position vectors. The local configuration defined by these molecules is described by the coefficients Cj+g,g G Q, i.e. 
by ( Cj+gi , Cj+g2 , . . . Cj+g^ , . . . Cj+g^ ) G (Z/n)'* , where /U is the number of molecules inside G- To each such multiplet 
of n'* vectors will correspond a relaxation time 

t( cj+gi , Cj+g2 > • • • Cj+g, , • • • q+g,. )• (17) 

Note that the quantities cj+g^ are only a small subset of the coordinates of c. Some of these relaxation times will 
be equal for symmetry reasons, but the notation remains generally valid. This definition holds at any site j, due 
to the translational invariance on C. The relaxation time Eq. (17) goes with jumps c c ± ej, i.e. it goes with 
-2^c,d + <5c,d+ej + <5c,d-ej , such that we have 

c^d c^+ej c^d_ej ^ ^^g^ 

jg£ '''( q+gi ' Cj+g2 ; ■ • • Cj+g< ; ■ ■ ■ Cj+gn ) 

If we further would like to introduce the possibility that molecules turn simultaneously, we will end up with 



p 



'''^^^ ( Q+gi ' Cj+g2 > • • • Cj+g« ) • • • Cj+g^ ) 

We can appreciate from such equations that the translational invariance on H is completely lost as the relaxation 
times depend on c. Although (atmost) "only" pn^ of them have been introduced, it is not so that we could define a 
unit cell of p sites on H and define a "Bravais lattice" of N/ {pn'^) cells in analogy with the matrix diagonalization 
procedure for a phonon problem. Indeed, there is translational invariance on C as all clusters j + G are congruent, 
but this does not carry over to H, e.g. ej+s ~ ej ^ ej+2 s — ^j+s, although (j+s)— j = (j-|-2s) — (j + s). 

The loss of translational symmetry has dramatic consequences. We can no longer diagonalize the jump matrix by 
simple arguments of translational symmetry and we run out of methods. There are other types of symmetry present, 
but these are not high enough to simplify the problem in any significant way. Therefore, in this case we cannot expect 
to recover a behaviour where the elastic intensity would vanish between the Bragg peaks, and the coherent quasielastic 
intensity would not reveal the correlations present. 
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